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Abstract 

We analyze the new threshold data for neutral pion photoproduction off pro- 
tons in the framework of heavy baryon chiral perturbation theory. We show 
that large loop corrections are needed to understand the S-wave multipole 
-Bo+ aiid that all pertinent low-energy constants can be understood within 
the framework of resonance exchange saturation. Previous inconsistencies in 
the description of this reaction in the threshold region are resolved. 
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Neutral pion photoproduction off protons has been a hot topic ever since the Saclay and 
Mainz 0] groups claimed a sizeable deviation from a so-called low-energy theorem (LET) 
for the electric dipole amphtude Eq+ derived in 1970 However, reexaminations of 
these data seemed to bring the empirical value in agreement with the theoretical prediction, 
see e.g. [|]. On the theoretical side, it was shown that the low-energy theorem of Refs. 
1^ is indeed incomplete and that the expansion of the electric dipole amplitude in 
powers of fi = Mt^/itl (with M^^ and m the pion and the nucleon mass, respectively) is 
slowly converging and therefore hard to pin down accurately. The numerical closeness of 
the empirical value of Eq^ at threshold with the one based on the incomplete LET has led 
to a flurry of proposals to reinterpret or resurrect the latter (for a detailed discussion, see 
e.g. [H). Two new developments, however, allow us to show in this letter that indeed there 
is no mystery about the threshold data for 7p — n^p if one performs a sufficiently accurate 
calculation in chiral perturbation theory. First, the theoretical framework to do just that 
was laid out in Ref. 0, as discussed very briefly below. Second, the new data from the 
TAPS collaboration have now been released |10[ and they show some discrepancies to the 



previously considered best data of Beck et al. |TT| . 

Let us briefly review the pertinent results of Ref. ^ . In that paper, heavy baryon chiral 
perturbation theory [|T3| was used to calculate the S-wave multipole i?o+ to order g"^ (where q 
denotes a small momentum) and the P-wave multipoles -Pi,2,3 to order [12|. The pertinent 
effective Lagrangian takes the form 

where the superscript (i) refers to the number of derivatives or meson mass insertions. The 
structure of Cj^n is discussed in detail in the review and a pedagogical introduction can 



be found in |jT5|. Besides the loop and pseudovector (pv) Born contributions, there are to 



this order two counter terms in the S-wave and one in the P-wave P3, 

Eo+ (io) = £;o^;™ (io) + (io) + ea,ujM^ + ea2Uj\ (2) 



P,{u) = P^°^-{u) + Pj''°^{u) 2 = 1,2, (3) 

Ps{u)=Pi°'''{u) + ebpu\q\ (4) 

with uj the pion energy in the cms system, q the pion momentum and e^/47r = 1/137.036. 
We have not made explicit the scale dependence of the low-energy constants ai and 02. 
In what follows, we use X = m. At threshold, we have Uq = M^o = 134.97 MeV. The pv 
Born contributions include the coupling proportional to the anomalous magnetic moment 
of the proton, Kp. In the chiral counting, these stem from the dimension two Lagrangian 
C^^N- Based on the data of Ref. 0, the three low-energy constants 01^2 and bp could be 
determined by a best fit. First, the numerical value for bp can be estimated from resonance 
exchange [1^, in this case from the A (in the static isobar model) and the vector mesons. 



V = p'' + uj, 

b'^'° = 6^ + 6^ = (9.7 + 3. 1) GeV"^ (5) 
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which shows that the vector meson contribution can not be neglected [0. The fitted value 
for bp is close to the number given in Eq.(^. However, letting the two S-wave constants ai 
and a2 completely free, they turn out to be very large in magnitude but of different sign. If 
one restricts these coefficients again by resonance exchange, one can only vary the A off-shell 
parameters within some bounds [0] and finds much smaller values for ai and 02 (typically 
a factor 20 smaller than in the free fit). This signals that there are either enormous higher 
order corrections or that the strong energy dependence of Eq+^lj) as suggested by the data 
of is incorrect. It is important to note, however, that the sum ai + 02 is roughly the same 
in both procedures. In effect, if no anomalously large coefficients appear, only this sum plays 
a role (in the threshold region). We also remark that the form of Eqs.(^ has lead to novel 
P-wave LETs for Pi and P2- These will be tested directly in polarization measurements at 
MAMI soon (for a somewhat model-dependent analysis, see |T9|). 



We can now use this formalism to analyze the new TAPS data of Fuchs et al. |T^. In 
Fig.l, we show the fit constrained by resonance exchange for the differential cross sections 
and in Fig. 2 for the total cross section. For the differential cross sections, all data up to 

= 160 MeV were used in the fit but only the ones up to E^ = 152 MeV are shown in 
Fig.l. Before discussing our fit parameters, we note that the differential cross sections above 
the 7i~^n threshold show much less of the pronounced bell shape as inferred from the older 
data and that the total cross section is somewhat decreased, which is of particular relevance 
for the extraction of Eq^. Let us now discuss in more detail the fits based on the theoretical 
framework of Ref. [^]. For the A, we keep the two 7A^A couplings fixed, gi = g2 = 5. 
The off — shell parameter Y is severely constrained by the A contribution to the magnetic 
polarizability of the proton, we set 0.1 < F < 0.14 , i.e. 6.4 < 6Pp < 7.5 (in units of 
10~^ fm^) [^. Also, Z is bounded by the A contribution to the vrA^ scattering volume 
033, —0.4 < Z < —0.2 X is varied within the range given in |TB|. We find X = 2.75 



23| , Y = 0.10 and Z = —0.21 which translates into 

(01 + 02)^"*"^ = (2.67 + 3.92) GeV-^ = 6.59 GeV-\ 651'° = 13.0 GeV^^ (6) 

The x^/dof is 2.21. We remark that the value for bp does indeed nicely agree with the 
resonance saturation estimate, Eq.(|^). Note that 6p'° does not depend on X and its possible 
values are strongly constrained by the ranges of Y and Z discussed above. It is thus gratifying 
that one can obtain such a consistent description of this low-energy constant. The sum 
(ai + 02) is consistent with the free fit value of 6.60 GeV~^ (we do not show the free fit since 
it is essentially the same as the resonance one). We also note that the resonance fit of Ref. 
0] already had ai + a2 = 6.67 GeV~^. The apparent mismatch between the free and the 
resonance fit discussed in has turned out to be an artefact related to the old data. 

In Fig. 3, we show the electric dipole amplitude £'0+. Its values at the n^p and the Tr+n 
threshold are 

Eo+(7r°j9) = -1.16 ■ 10-yM^+ , Eo+{n+n) = -0.44 • 10-VM^4- , (7) 



to be compared with Eof^n^p) = -1.31 ± 0.08 ■ 10-^/A/U+ and Eof{7r+n) ~ -0.4 



10~'^/M^+ (as read off from Fig. 4 of |T^). The value of ai + 02 in Eq.(|^) amounts to an 
i?o+-contribution of +0.3 from vector mesons and +0.4 from the A (in units of 10~^/M.,^+). 
Almost the same number for the sum of vector meson and nucleon resonance countributions 
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to £"0+ at threshold is reported in Ref. p4[ (see also |2^). In the threshold region, the shape 



for Eq^{uj) shown in Fig. 3 can be well represented by a two-parameter fit of the form (as 
discussed in some detail in Ref. 0) 

Eo+(^) = -a-6^1-^Va;2, (8) 

with uJc = 140.11 MeV the pion energy at the n^n threshold. We find a = 0.44 ■ 10~'^/Mt^+, 
6 = 2.9 ■ 10-VM^+ and bp = 12.9 GeV'^. The xV^of is 2.22, i.e. almost identical to the 
one of the resonance fit. This indicates that the mild slope of Re Eoj^{uj) behind the Tr+n 
threshold is not significant. The value for b is somewhat below the one estimated from the 
Fermi-Watson theorem, b^^ = 3.7 • 10^^ /Mt^j^^. Note, however, that this is based on the 
assumption of exact isospin symmetry, whereas the clearly visible cusp effect in Eq^{u) is 
due to the pion mass difference, i.e. an isospin-violating effect. A more consistent treatment 
of such effects is certainly needed. For a study of the cusp effect in Eq^ in terms of a multi- 
channel S-matrix, see Ref. ||2^. We note that the small value of a is a clear indication 
of chiral loops - with simple pseudovector Born terms one can not get such a small value 



(this was already stressed in 0). 7P — ir^p has also been remeasured at Saskatoon p7 
Between n^p and 7r~^n thresholds, the new SAL data are consistent with the ones of Ref. 
10| . For larger energies, however, the new SAL data agree with the older Mainz data 0. 



This does not affect the threshold value of Eq^ but rather leads to a larger value of 6p. The 
experimental discrepancy remains to be clarified. 

For the respective slopes of the P-wave multipoles, the LETs together with the best res- 
onance fit give Pi/\q \ = 0.480 GeV-^, P^/lql = -0.512 GeV-^ and Ps/\q \ = 0.544 GeV-^, 
i.e. all P-waves -Pi,2,3 are of the same magnitude close to threshold. Consequently, the 
photon asymmetry S(^) to be measured at MAMI is expected to be small in the threshold 
region since S ~ (I-P3P — l-P2p)- 

To summarize, we have shown that within the framework of chiral perturbation theory, 
one is able to consistently understand the new threshold data for the reaction 7p — vr^p. 
Loop effects are clearly visible in the S-wave. Furthermore, the three low-energy constants 
are fully understood within the framework of resonance saturation. In contrast to common 
folklore, these resonances pose no problem and do not have to be treated as dynamical 
degrees of freedom (as long as one stays in the threshold region). In the next step, this 
formalism should be extended to electroproduction to discuss the new data from NIKHEF 
and MAMI m. 
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FIGURES 



Fig.l Differential cross sections in the threshold region (in nb/sr) for lowest 9 values of the 
photon lab energy versus the cm scattering angle 9. The solid line is the best 



resonance fit, the data are from |TD 



Fig. 2 Total cross section in the threshold region (in /^b) versus E^. For notations, see Fig.l. 
Fig. 3 The real part of the electric dipole amplitude in the threshold region. 
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